The renormalization group and quark number fluctuations in the Polyakov loop 
extended quark-meson model at finite baryon density 
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Thermodynamics and the phase structure of the Polyakov loop-extended two flavors chiral quark- 
meson (PQM) model is explored beyond the mean-field approximation. The analysis of the PQM 
model is based on the functional renormalization group (FRG) method. We formulate and solve 
the renormalization group flow equation for the scale-dependent thermodynamic potential in the 
presence of the gluonic background field at finite temperature and density. We determine the phase 
diagram of the PQM model in the FRG approach and discuss its modification in comparison with 
the one obtained under the mean-field approximation. We focus on properties of the net-quark 
number density fluctuations as well as their higher moments and discuss the influence of non- 
perturbative effects on their properties near the chiral crossover transition. We show, that with an 
increasing net-quark number density the higher order moments exhibit a peculiar structure near the 
phase transition. We also consider ratios of different moments of the net-quark number density and 
discuss their role as probes of deconfinement and chiral phase transitions. 
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I. INTRODUCTION 

Thermodynamic properties of a strongly interacting 
matter at nonzero baryon density and high temperature 
were quantified numerically within Lattice QCD (LQCD) 
. The LQCD results demonstrate that QCD exhibits 
both dynamical chiral symmetry breaking and confine- 
ment at finite temperature and densities. The LQCD 
equation of state indicates a clear separation between 
the confined hadronic phase and deconfined phase of 
quark-gluon plasma. However, since the thermodynam- 
ics at large baryon densities and near the chiral limit 
is still not accessible from the first principle LQCD cal- 
culations, many phcnomcnological models and effective 
theories have been developed |4H18j. 

The hadronic properties at low energy as well as the 
nature of the chiral phase transition at finite tempera- 
ture and densities have been successfully explored and 
described in such effective models. The physics of color 
deconfinement and its relation to the chiral symmetry 
breaking has been recently studied in terms of effec- 
tive models. The idea to extend the existing chiral La- 
grangians such as the Nambu-Jona-Lasinio or the quark- 
meson, by introducing coupling of quarks to uniform tem- 
poral background gauge fields (Polyakov loops) was an 
important step forward in these studies 0, [IH ■ 

It was shown that the Polyakov loop extended Nambu- 
Jona-Lasinio (PNJL) @ or quark-meson (PQM) 
models can reproduce essential properties of the QCD 
thermodynamics obtained in the LQCD already within 
the mean-field approximation. However, to correctly 
account for the critical behavior and scaling properties 



near the chiral phase transition one needs to go beyond 
the mean-field approximation and include quantum fluc- 
tuations and non-perturbative dynamics. This can be 
achieved by using the methods based on the functional 
renormalization group (FRG) fT~9l - [25j ] . 
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Following our previous work [26| we use a truncation 
of the PQM model which is suitable for the functional 
renormalization group analysis to describe the thermo- 
dynamics beyond the mean-field approximation. The 
functional renormalization group approach in the PQM 
model is used to take into account fluctuations of the me- 
son fields, while the Polyakov loop is treated as a back- 
ground field on the mean-field level. In contrast to the 
previous work [26j we extend our calculations to the fi- 
nite chemical potential. We determine the phase diagram 
and the position of the critical end point (CEP) in the 
PQM model by exploring the dependencies of the chi- 
ral order parameter and the quark number susceptibility 
on thermal variables. We calculate the moments (cu- 
mulants) of the net-quark number density fluctuations 
(c„) at finite temperature and chemical potential in the 
presence of mesonic fluctuations. We discuss the influ- 
ence of non-perturbative effects on properties of the first 
four c n -moments near the chiral crossover transition. We 
show that c n -cumulants exhibit a peculiar structure and 
for sufficiently large values of the chemical potential can 
be negative near to the crossover transition. We calcu- 
late the ratios Czjc\ and c^jci and discuss their roles as 
probes of deconfinement and chiral phase transitions. We 
summarized properties of different susceptibilities near 
the chiral phase transition at finite net-quark density 
within the Landau mean-field and the scaling theories. 
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II. THE POLYAKOV-QUARK-MESON MODEL 

The model which is used in this paper to explore the 
chiral phase transition at finite temperature and density 
within the FRG approach is the Polyakov loop-extended 
two flavor quark-meson model. In general, the PQM 
model, being an effective realisation of the low-energy 
sector of the QCD, cannot describe confinement phenom- 
ena because the local SU (N c ) invariance of the QCD is 
replaced by the global symmetry. However, it was argued 
that by connecting the chiral quark-meson (QM) model 
with the Polyakov loop potential the confining properties 
of QCD can be approximately accounted for J7l [H|, EBI- 

The Lagrangian of the PQM model reads [15| 



<1 WD, 



1 



g(a + i^fif)] q + -(c^cr) + -(9 M 7r) 



-C/((T,7?) -U{i,t) . 



(1) 



The coupling between the effective gluon field and quarks 
is implemented through the covariant derivative 



Dn — dn — iAu, 



(2) 



where = gA a ^\ a /2. The spatial components of the 
gluon field are neglected, i.e. A^ = S^qA . Moreover, 
U(£,£*) is the effective potential for the gluon field ex- 
pressed in terms of the thermal expectation values of the 
color trace of the Polyakov loop and its conjugate 



f 



<Tr c L(2)>, t 



1 



(Tr c Lt(a?)>, (3) 



with 



L(x) = 'Pexp 



dr A4 (2?, r) 



(4) 



where V stands for the path ordering, /3 = 1/T and 
A4 = iAq. In the 0(4) representation the meson 
field is introduced as <fr = (a, tt) and the correspond- 
ing SU(2)l SU(2)r chiral representation is defined by 
a + if ■ 7?75 . 



The purely mesonic potential of the model, U (a, 7?), is 
defined as 



[/(a^)^(a 2 



-*2 2 \ ^ 

7T — V ) — C(T, 



(5) 



while the effective potential of the gluon field is 
parametrized in such a way as to preserve the Z{3) in- 
variance: 



U{IJ*) b 2 (T) 
T 4 2 



The parameters, 



+ t z ) + ^{tt) 2 . (6) 



a + ai 



(7) 



with a = 6.75, ai = -1.95, a 2 = 2.625, a 3 = -7.44, 
63 = 0.75 and 64 = 7.5 were chosen to reproduce the 
equation of state of the pure SU C (3) lattice gauge theory. 
Consequently, at To ~ 270 McV the potential © yields 
the first-order deconfincment phase transition. Several 
alternative parametrizations of the effective potential of 
the gluon field were also explored in Refs. |27H29| . 



A. The FRG method in the PQM model 

In order to formulate a non-perturbative thermody- 
namics within the PQM model we adopt a method based 
on the functional renormalization group (FRG). The 
FRG is based on an infrared regularization with the mo- 
mentum scale parameter of the full propagator which 
turns the corresponding effective action into the scale 
fc-dependent functional Tk (T9T - [22l |. 

In the PQM model the formulation of the FRG flow 
equation would require implementation of the Polyakov 
loop as a dynamical field. However, in the current cal- 
culation we treat the Polyakov loop as a background 
field which is introduced selfconsistcntly on the mean- 
field level. 



Following our previous work [26j we formulate the flow equation for the scale-dependent grand canonical potential 
for the quark and mesonic subsystems 



12tt 2 }E n 

4N c Nf 
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Here ns(£' 7rj(T ; T) is the bosonic distribution function 



cxp(£™/T) - 1 
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FIG. 1: The phase diagrams for the PQM model in the mean-field approximation (left panel) and in the FRG approach (right 
panel). The shaded regions are defined by 5%-deviations of the temperature derivative of the chiral order parameter, da/dT, 
from its maximal value. The arrows show the lines corresponding to different values of j-i/T. The dashed curves indicate 
isentropes for s/n q — 2, 5, 10, 25, 50, 200. 



with the pion and sigma energies given by 

K = \lk 2 +n' k , E a 



k2 + n k + 2 P n k ; 



where the primes denote derivatives with respect to p and Q, = + ca. The functions N(£, £*; T, /j.) and N(£, £*; T, fi) 
defined by 



N(£ t . T u) = 1 + 2t exp[P(E q -»)]+£ expjg^gg - rf\ 

[ ' ' ' W 1 + 3* exp[2/3(E q - p)} + 3£* exp[(3{E q - fj,)} + exp{3f3{E q - p:)} ' 
N(£,£*;T,f,) = N(t,£;T,-») 

are fermionic distributions modified because of the coupling to gluons. The quark energy is defined by 



E q = Vfc 2 + 2 5 2 (0 - 



(9) 
(10) 

(11) 



r 



The minimum of the thermodynamic potential is de- 
termined by the stationarity condition 



da 



dfl k 
da 



(12) 



il k ^o(£, £*; T, fi), as the function of the Polyakov loop 
variables £ and £*. The full thermodynamic potential 
Vl(£,£*;T, fi) in the PQM model which includes quarks, 
mesons and gluons degrees of freedom is obtained by 
adding to fl k ->o(£, £*; T, fi) the effective gluon potential 
U(£,£*): 



The flow equation ([8]) is solved numerically with the 
initial cutoff A = 1.2 GeV (see details in Ref. HI). The 
initial conditions for the flow are chosen to reproduce 
the vacuum properties: the physical pion mass — 
138 MeV, the pion decay constant = 93 MeV, the 
sigma mass m a = 600 MeV, and the constituent quark 
mass m q = 300 MeV at the scale k = 0. The symmetry 
breaking term, c = m^f^, corresponds to an external 
field and consequently does not flow. In this work we 
neglect the flow of the Yukawa coupling g, which is not 
expected to be significant for the present studies (see e.g. 
Rcfs. 0,|3l|). 

By solving the equation ([S]) one gets thermody- 
namic potential for quarks and mesonic subsystems, 



Sl(e, t;T, n) = n k ^ (£, £*■ T, /x) + U{i, £*). 



(13) 



At a given temperature and chemical potential, the 
Polyakov loop variables, £ and £*, are determined by the 
stationarity conditions: 

d -n(£,£*;T, f i)=0, (14) 



d£ 
d_ 

at 



n(£,£*;T,») = 0. 



(15) 



The thermodynamic potential (|13[) does not contain 
contributions of statistical modes with momenta larger 
than the cutoff A. In order to obtain the correct high- 
temperature behavior of thermodynamic functions we 
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FIG. 2: The baryon number density normalized by T 3 , ci = n q /T 3 , as a function of temperature for different values of fi/T 
for the PQM model in the mean-field approximation (left panel) and in the FRG approach (right panel). 



need to supplement the FRG potential with the contri- 
bution of high-momentum states. A simplified model for 
implementing such states was proposed in Ref. [HJ where 
the contributions of the k > A states to the flow is ap- 
proximated by that of a non- interacting gas of quarks and 
gluons. For the PQM model we generalize this procedure 
by including for k > A the flow equation of interacting 
quarks with Polyakov loops [26[ 



N c N f k 3 



3tt 2 



(16) 



l-N(£,£*;T,fi)-N(£,t;T,fi) , 



where the dynamical quark mass is neglected. In addi- 
tion, since the effective gluon potential U (£,£*) is fitted 
to reproduce the Stefan-Boltzmann limit at high tem- 
peratures, the explicit gluon contribution is omitted for 
consistency. 



J 



To obtain the complete thermodynamic potential of 
the PQM model we integrate Eq. (|T6| from k = oo to 
k = A where we switch to the PQM flow equation ©. Di- 
vergent terms in the high-momentum flow equation (fTf))) 
are independent of mesonic and gluonic fields as well as of 
temperature and chemical potential. Consequently, such 
terms can be absorbed to an unobservable constant shift 
of the thermodynamic potential. 



B. The mean-field approximation 

To illustrate the importance of mesonic fluctuations 
on thermodynamics in the PQM model we will compare 
the FRG results with those obtained under the mean- 
field approximation. In the latter, both quantum and 
thermal fluctuations are neglected and the mesonic fields 
are replaced by their classical expectation values. 



In the PQM model the thermodynamical potential derived under the mean-field approximation reads [if 

fl MF = U(£,£*) + U((a),(ir) = 0) + n gg -((a), £,£*). 
Here, the contribution of quarks with dynamical mass m q = g{a) is given by 



where 



n qS ((a),£,t 



-2N f T 



d 3 p ( N C E, 



(27T) 



T 



q - + In g (+) «<r>, £, t ; T, M ) + In <?(") ( (a) ,£,£*; T, ft) 



(17) 
(18) 

(19) 
(20) 

and Eg = ^Jp 2 + is the quark quasi-particle energy. The first term in Eq. (|18[) is a divergent vacuum fluctuation 

contribution which has to be properly regularized. Following our previous studies (33j we use the dimensional regu- 
larization, which introduces an arbitrary scale M. The renormalization procedure allows to compensate the physics 
dependence on this scale. The finite contribution of the vacuum term to Eq. (fTS)) reads [33[ 

N c Nf 



g^ + >{(a),e,e*;T,fi) = 1 + 3£cxp[-(E q - ft)/T] + U* exp[-2(E q - fi)/T] + exp[-3(E q - fi)/T], 
g(-\{a),£,t;T,fi) = g^((a),£*,£;T,-ft); 



r\vac 



8tt 2 



(21) 
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FIG. 3: The coefficient C2 as a function of temperature for different values of fi/T for the PQM model in the mean-field 
approximation (left panel) and in the FRG approach (right panel). 



The importance and influence of this contribution on the thermodynamics of chiral models was demonstrated and 
studied in detail in Refs. [H[ and [35j . 



The equations of motion for the mean fields are ob- 
tained by requiring that the thermodynamic potential is 
stationary with respect to changes of (a), £ and £*: 

9Qmf _ QVLmf _ 9Q.MF _ „ 

The model parameters are fixed to reproduce the same 
vacuum physics as in the FRG calculation. 



III. THERMODYNAMICS OF THE PQM 
MODEL 

The thermodynamic potential introduced in Eqs. (fl~3]) 
and (fT7|) provides basis to explore critical properties of 
the PQM model at finite baryon density within the FRG 
approach and under the mean-field approximation. 

To find the potential at finite temperature and chem- 
ical potential one needs to solve the FRG flow equation 
flSJ). This equation is solved by numerical methods based 
on the Taylor series expansion [2(| |34[ . This method is 
successful in stud ying thermodynamics at finite density 
and temperature |24l |2(I [35j in the regime where the 
system exhibits the crossover or the second-order phase 
transition. For the solution of the FRG flow equations 
in the regime of the first-order phase transition, where 
the thermodynamical potential develops two minimums, 
different numerical methods are required [H, [36|. In 
the present work we restrict our considerations only to 
the parameter range where the PQM model exhibits the 
crossover or the second-order chiral phase transitions. 



A. The phase diagram 

The PQM model is expected to belong to the same uni- 
versality class as QCD and thus should exhibit a generic 
phase diagram with the critical point at non-vanishing 
chemical potential. In the chiral limit the phase diagram 
is identified by divergences of the chiral susceptibility. At 
finite quark mass the chiral transition is of crossover type. 
In this case the pseudocritical temperature and chemi- 
cal potential are located by determining the maximum of 
the chiral susceptibility or the temperature derivative of 
the chiral order parameter. The position of the CEP is 
identified by the properties of the sigma mass. The tem- 
perature and chemical potential where the sigma mass 
vanishes correspond to the position of the CEP. One can 
equivalently consider the net-quark number fluctuations 
to identify the critical point which according to Z(2) uni- 
versality diverge at the CEP. 

Fig. [T] shows the phase diagrams of the PQM model 
obtained in the FRG approach and in the mean-field ap- 
proximation. For the physical pion mass and moderate 
values of the chemical potential the PQM model exhibits 
a smooth crossover chiral transition. In Fig. [T] we de- 
fine the transition region as bands where the temperature 
derivative of the order parameter exhibits 5%-dcviations 
from its maximal value. At larger chemical potentials the 
crossover line terminates at the CEP where the transition 
is of the second order and belongs to the universality class 
of the three dimensional Ising model. 

Comparing the resulting phase diagrams of the PQM 
model obtained within the mean-field approximation and 
FRG approach we find a clear shift of the position of 
the chiral phase boundary to higher temperatures due to 
mesonic fluctuations. Such shifts were previously found 
in the QM model within the FRG approach |23|, l35j . 
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FIG. 4: The coefficient C3 as a function of temperature for different values of n/T for the PQM model in the mean-field 
approximation (left panel) and in the FRG approach (right panel). 



However, in our studies due to the gluonic background, 
which is explicitly included in our FRG calculations, we 
also find a significant shift of the CEP to higher temper- 
ature. 

In the following we focus on observables that are re- 
lated to the non- vanishing net-quark number density and 
consider moments of quark number fluctuations. We dis- 
cuss how such fluctuations depend on the quark chemi- 
cal potential in the presence of the gluonic background 
within the FRG approach. 



B. Quark number density fluctuations 

The fluctuations of conserved charges are observables 
that provide direct information on critical properties re- 
lated with the chiral symmetry restoration. Fluctuations 
related with the baryon number conservation are of par- 
ticular interest. In an equilibrium medium a divergence 
of the net-quark number susceptibility is a direct evi- 
dence for the existence of the CEP J40|. Consequently, 
any non-monotonic dependence of these fluctuations on 
collision energy in heavy-ion collisions was proposed as a 
phenomenological method to verify the CEP [40j . In a 
non-equilibrium system, the net-quark number suscepti- 
bility was also shown to signal the first-order chiral phase 
transition due to spinodal decomposition [Io| . 

The fluctuations of the net-quark number density are 
characterized by the generalized susceptibilities, 



c„(T) = 



g"[p(?»/r 4 

0(M/r)» 



(23) 



The first coefficient c\ = n q /T 3 characterises the re- 
sponse of the net-quark number density to changes in 
the quark chemical potential. The second coefficient 



C2 



Xq_ 

J^2 



1 



VT 3 



(24) 



with SN q = N q — (N q ) is proportional to the susceptibil- 
ity of the net-quark number density, \ q . The third and 
fourth order moments can be expressed through SN q as 



1 



f"3 



f'4 



VT 3 
1 

VT 3 



((SNq) 3 ), (25) 

({(SNq^-SdSNqff). (26) 



The coefficients c n (T) are sensitive probes of the chiral 
phase transition. They indicate the position, the order 
and, in case of the second-order phase transition, the 
universality class of corresponding phase transition. The 
net-quark number density, n qi is discontinuous at the 
first-order transition, whereas the susceptibility C2 and 
higher cumulants diverge at the CEP (37H39j . In the chi- 
ral limit and at non-zero chemical potential, all general- 
ized susceptibilities c n (T) with n > 2 diverge at the 0(4) 
chiral critical line [4l[. Moreover, they also diver ge a t the 
spinodal lines of the first-order phase transition |K| . 

A very particular role is attributed to the so-called kur- 
tosis of the net-quark number fluctuations j4lTj43l | which 
is defined as the ratio 



R4.2 — — 



(27) 



This key observable is not only sensitive to the chiral 
but also to the deconfinemcnt transition. At vanish- 
ing chemical potential, in the asymptotic regime of high 
and low temperatures the kurtosis reflects quark content 
of the baryon-number carrying effective degrees of free- 
dom [4l|, |43|. Therefore, at low temperatures in the con- 
fined phase, i?4 ; 2 — Nq = 9 while for high temperatures 
one recovers an ideal gas of quarks with R^i ~ 1 

The properties of different moments of the net-quark 
number fluctuations in the presence of the chiral phase 
transition were studied in the literature both in terms of 
the LQCD [1| as well as in different models [U H, H, 
HH, 13 EH- In particular, the importance of the quark 
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FIG. 5: The coefficient C4 as a function of temperature for different values of /i/T for the PQM model in the mean-field 
approximation (left panel) and in the FRG approach (right panel). 
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FIG. 6: The kurtosis i?4,2 as a function of temperature for different fi/T calculated in the PQM model under the mean-field 
approximation (left panel) and in the FRG approach (right panel). 



susceptibility and the kurtosis as signatures of the decon- 
finement and the chiral phase transition as well as the 
CEP was discussed [43}. The influence and dependence 
of these fluctuations on the quark mass was also analyzed 
on the lattice and in effective models [H [H, [H, [H 53 . 
However, only little is known about chemical potential 
dependence of the higher cumulants c„ , particularly with 
n > 2. Such dependence can be obtained in the PQM 
model from the thermodynamic potential introduced in 
Eqs. (fT"3")) and (fTTj) . In Figs.[2HS]we quantify the first four 
moments obtained in the PQM model under the mean- 
field approximation and in the FRG approach for differ- 
ent values of the ratio /i/T. The lines of the fixed ratio 
/i/T also indicated on the phase diagram in Fig. [TJ 

Fig. [2] shows the net-quark number density normalized 
by T 3 , ci = n q /T 3 , for different values of /i/T. The cu- 
mulant c\ as well as all C2 n +i, for n = 1, 2, 3, . . ., are van- 
ishing at zero chemical potential, /i = 0. At finite \x and 
in the chiraly broken phase the net-quark number, n q , is 
strongly increasing function of fi/T. In the low tempera- 
ture phase due to the large dynamical quark mass the n q 



increases approximately as sinh(3/i/T). This is a direct 
consequence of the "statistical confinement" properties 
of the PQM model which for small expectation values 
of Polyakov loops I 1 suppresses the single and dou- 
ble quark modes in the partition function as seen from 
Eqs. © and (fTT)|). At high temperatures/chemical poten- 
tials n q behaves as polynomial in fx/T. For a fixed \xjT 
and in the pseudo-critical region where the chiral sym- 
metry is partially restored there is a clear change in the 
temperature dependence of n q . At the pseudo-critical 
temperature the density exhibits a kink-like structure 
which is particularly evident in the mean-field calcula- 
tions at larger n/T. In the FRG calculations and at the 
corresponding /i/T the above kink-like structure in the 
density n q is suppressed. Thus, the meson fluctuations 
quantified by the FRG method provide smoothing of the 
net-quark density near the crossover transition. 

The appearance of the crossover chiral transition in 
the PQM model is also transparent when considering the 
structure of isentropic trajectories. Because of the en- 
tropy and baryon-numbcr conservation such trajectories 



8 




15 



T 5 



150 200 250 
T [MeV] 



300 



1 . ' 


i 1 i 1 


- s/n b =200 ' 

s/n b =50 

s/n b =25 

" s/n b =10 

■ s/n b =5 

i 




i.i. 



100 150 200 250 
T [MeV] 



300 



FIG. 7: The temperature dependence of the kurtosis calculated in the PQM model at fixed values of the entropy density 
to net-quark density (s/n q ) under the mean-field approximation (left panel) and in the FRG approach (right panel). 



correspond to contours of the constant entropy density 
per net-quark number density, s/n q , in the temperature- 
chemical potential plane. Fig. Q] shows isentropes for 
different s/n q in the PQM model obtained under the 
mean-field approximation and in the FRG calculations. 
There is a clear change of slopes of isentropes along the 
transition line indicating the change of the equation of 
state. The qualitative behavior of contours of constant 
s/n q obtained in the PQM model is similar to that calcu- 
lated previously within the QM model [35|. In particular, 
the isentropes remain smooth when the effect of long- 
wavelength meson fluctuations is consistently included 
in the presence of the gluon background. There are also 
no indications of any focusing towards the CEP in the 
PQM model. 

The influence of the finite chemical potential on c 2 
(which is proportional to the net-quark number sus- 
ceptibility C2 = Xq/T 2 ) is shown in Fig. [3] for the 
mean-field approximation and the FRG approach. At 
vanishing chemical potential the cumulant 02 increases 
monotonously with temperature. However, at finite 
chemical potential, the susceptibility C2 develops a peak 
structure. The amplitude of this peak increases with the 
chemical potential towards the CEP where C2 diverges. 
In the high temperature/chemical potential phase C2 con- 
verges to the Stefan Boltzmann limit 

(28) 

As it is seen in Fig. [3] the peak structure in C2 is more 
pronounced in the mean- field approximation at fi/T = 1 
than in the FRG at /i/T = 1.5. This is in spite of the 
fact that the location of the CEP in the FRG is closer 
to the line of [i/T = 1.5 than the corresponding one for 
the mean-field approximation to /i/T = 1 (sec Fig. [T|). 
This shows that the criticality of C2 as a function of the 
distance to the CEP appears earlier in the mean-field 
approximation than in the FRG approach. This result is 
in agreement with the previous studies in the QM model 



showing that the critical region shrinks due to mcsonic 
fluctuations p3| . 

The cumulants c\ and ci arc sensitive to changes in 
the chemical potential and are influenced by the meson 
fluctuations and the gluon background. However, both 
these coefficients are positive for all values of /i and T. 
At finite chemical potential the positivity of the first two 
coefficients is not preserved for c n -moments with n > 2. 
Figs. 0] and O show the third and the fourth order cumu- 
lants of the net-quark number density for different values 
of /i/T. For a large /i/T both these susceptibilities ex- 
hibit a peculiar structure in the transition region. There 
is a minimum of C3 which can reach negative value. The 
amplitudes of this minimum is strongly suppressed by 
meson fluctuations in FRG approach. 

The fourth order cumulant is strictly positive for van- 
ishing chemical potential. However, for higher values of 
jU/T, C4 becomes negative in the vicinity of the crossover 
transition. This manifests the broadening of the net- 
quark number distribution in comparison to the Gaussian 
one. Large values of C4 in the broken phase infer that the 
distribution becomes narrower than the Gaussian. The 
chemical potential independent Stcfan-Boltzmann limit 
cf B = 2N c Nf/n 2 is reproduced at temperatures T T c . 

Comparing the mean-field with the FRG results for 
both C3 and C4 we conclude that mesonic fluctuations es- 
sentially modify properties of different generalized quark 
susceptibilities. In the transition region C3 and C4 are 
suppressed in the FRG relative to the mean-field results. 
Thus, the mean-field approach is only an approximate 
method to describe static thermodynamics near the chi- 
ral phase transition. 

We have already indicated that the ratios of different 
c n are sensitive to the deconfincmcnt and the chiral phase 
transitions. Figs. [5] and [7] show the kurtosis R^2 = C4/C2 
calculated as a function of temperature along different 
paths in the temperature-chemical potential plane quan- 
tified by fixed ji/T and sjn q . In the PQM model the 
kurtosis drops from R^2 — 9 to R^2 ~ 1 in the transi- 
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FIG. 8: The ratio of C3 to c\ as a function of temperature for different values of fi/T obtained in the PQM model under the 
mean-field approximation (left panel) and in the FRG approach (right panel). 



tion region. Such a change is explicitly attributed to the 
change in quark content of baryon number carrying effec- 
tive degrees of freedom [4l|. In the PQM model and in 
the chiraly broken phase the effective three quark states 
dominate, while at high temperatures single quarks pre- 
vail. In the mean-field approximation the kurtosis ex- 
hibits a peak at the transition temperature at fj, = 0. 
The height of these peaks depends not only on the pion 



26, 33, 43, 47] but also on the value of the chemi- 



cal potential. The mesonic fluctuations weaken the peak 
structure both at finite and at vanishing quark density. 
For finite chemical potential the kurtosis becomes nega- 
tive following the same trends as seen in the fourth order 
cumulant. 

A direct information on quark content of baryon car- 
rying effective degrees of freedom in the low temperature 
phase is also contained in the c^/ci ratio. At low tem- 
perature where thermodynamics is dominated by three- 
quark modes c^/ci = 9 for any value of chemical po- 
tential. At low temperatures and at zero chemical po- 
tential ci and C3 vanishes but their ratio is finite. At 
asymptotically large temperatures and for /i — y this 
ratio diverges. The ratio c^/c\ similarly as the kurto- 
sis -R4.2 = C4/C2 exhibits strong variations in the phase 
transition region. It develops a peak with height which 
increases with /i/T and for sufficiently large /i/T it devel- 
ops a deep structure. This is the case for both mean-field 
and FRG calculations, however variations at correspond- 
ing /i/T in the FRG are strongly suppressed because of 
meson fluctuations. 



C. Scaling properties of generalized susceptibilities 
at finite chemical potential 

The general trends and behavior of different suscepti- 
bilities calculated in the last section can be understood 
considering their scaling properties near the chiral tran- 
sition. Under the mean-field approximation such scaling 



can be inferred from the Landau theory of phase tran- 
sitions, where the singular part of the thermodynamic 
potential is a polynomial in an order parameter a, 



'(7>;a) 



ha. 



(29) 



with h being a symmetry breaking term. 

For the chemical potential much smaller than the crit- 
ical temperature T c (/i = 0) of the second-order chiral 
phase transition the mass term a(T, /z) is parameterized 
as 



a(T,ri=A-(T-T c ) + B 2 n 2 



(30) 



where both coefficients A and B2 are positive. In general, 
the effective quartic coupling b > 0, is also T and /i- 
dependent. However, this dependence is irrelevant near 
the critical line a(T,fi) = and away from the CEP or 
the tri-critical point (TCP). Therefore in this case we can 
also neglect the higher order polynomial contribution by 
setting c = 0. In the chiral limit h = 0, the critical 
properties of c-i and C4 are obtained from Eq. (|29p : 



(T-T C )9(T C -T), 



AB 2 
bT 2 

^9{T C -T). 



(31) 
(32) 



Thus, C2 is not diffcrentiable at the critical temperature, 
while C4 has a discontinuity. 

For a finite quark mass i.e. for h ^ 0, the second- 
order transition is turned to the crossover and the sharp 
structures in C2 and C4 are smoothened. Therefore in 
the PQM model, the peak structure appearing in C4 is 
directly linked to the quark mass and closely connected 
to the partial restoration of the chiral symmetry. From 
Eq. (|32[) it is clear, that the kurtosis i?4,2, driven by the 
C4 has a discontinuity at T c . 
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In the FRG approach the critical behavior of gener- 
alized susceptibilities obtained under the mean-field ap- 
proximation will be modified by the long wave meson 
fluctuations. Detailed studies in the QM model showed 
that the FRG method can correctly account for long- 
range correlations resulting in the 0(4) critical behavior 
of thermodynamic functions (24| . The gluon background 
is not modifying critical dynamics related with the chiral 
symmetry. Therefore in the chiral limit the PQM model 
belongs to the 0(4) universality class. The singular part 
of the thermodynamic potential is controlled by the crit- 
ical exponents of the three-dimensional 0(4)-symmetric 
spin system. At vanishing chemical potential 



n smg 



(T-T c 



(33) 



leading to the following scaling of generalized suscepti- 
bilities 



C2n+1 



(T - T c ) 




2—n—a 



(34) 
(35) 



for n — 1,2,3,... 

In the mean-field approach the critical exponent a = 0. 
The quantum fluctuations within FRG renormalize the 
exponent to a ~ —0.21 expected in the 0(4) universality 
class. Therefore fluctuations lead to weakening of singu- 
larities. The finite quark mass further smooths the tem- 
perature dependence of C2, C3 and C4 as seen in Figs. |3l |4] 
and [5] The kurtosis follows singular behavior of C4 and 
for /i = exhibits a step- like structure at T c . 

For finite chemical potential, but still away from the 
CEP or the tricritical point (TCP), the coefficient a(/x, T) 
in the Landau potential can be parameterized as 



a(T, fi) = A ■ (T - T c ) + B x ■ (/x - /x c ) 



(36) 



while the quartic coupling b > and we still keep c = 0. 

In this case one gets the following expressions for sus- 
ceptibilities 



^sing 



^sin g 



B x a 
2T 3 b 

2bT 2 

sing 



6(-a), 
t-a), 



0. 



(37) 

(38) 
(39) 



Thus, C2 exhibits a discontinuity while the singular part 
of C4 is vanishing along with C3 and higher order cumu- 
lants. However, there is the next to leading order con- 
tribution to a(T, n) owing to non- vanishing curvature of 
the transition line in the (T, /x)-plane. In order to take 
it into account one needs to add an extra contribution 
B2 ■ (/i — /x c ) 2 to the right-hand side of Eq. (|36|). In this 
case the c™ 9 behaves as in Eq. (f3"2")l while the leading 
contribution to c™ 9 and c^™ 9 have the structure as in 
Eq. (|37| but with a modified a(T, [i). The third order cu- 
mulant also develops nonzero values in the broken phase 



3-B2-B1 

bT 



kurtosis has a step-like behavior but here the curvature 
of the phase diagram controls whether the kurtosis is an 
increasing or decreasing function at the phase transition 
because of two competing step structures in c™ 9 and 

sinq 

Including quantum fluctuations as in the FRG ap- 
proach the above MF scaling is modified to the following 
0(4) relations 



n sms 

sing 



\2-ci 



(-a) 5 



(40) 
(41) 



resulting in divergence of the C3, C4 and all higher order 
cumulants along the 0(4) critical line with the specific 
heat critical exponent a ~ —021. The kurtosis is clearly 
divergent at T c = T((i c ) following the singularity of the 



„stng 



cumulant. 



The above scaling properties of the net-quark num- 
ber fluctuations at finite chemical potential are modified 
when approaching the TCP in the chiral limit ar CEP at 
finite quark mass. 

Close to the TCP the effective Landau potential has 
the structure as in Eq. (|29[) where the parameter a(T,fi) 
has a linear dependence on the reduced temperature and 
chemical potential 

a(7» = A a ■ (T - T c ) + B a ■ (/i — /i c ). (42) 

The quartic coupling tends to zero as 

b{T, n) = A b ■ (T - T c ) + B b ■ {n - Mc ). (43) 

The six-order coupling c > 0, as it is required by the 
stability of the theory. Consequently, within Landaus 
theory and with the above parametrization of the poten- 
tial coefficients one gets the following relations for the 
leading contribution to generalized susceptibilities, 



^sing 



B a -a 



|) (4c)™- 2 B™T n - 4 



2r(i) 



(b 2 - 4ac)"- 3 / 2 



(44) 

n > 1. (45) 



resulting in the divergent kurtosis at the TCP as: 



H 4,: 



(6 2 -4ac)" 



6(—a). As in the case of /i = 0, the 



From Eq. (|45|) one concludes that the c™ 9 is in- 
versely proportional to a distance from the TCP along 
the a(T, fx) = line. When approaching the TCP from 
any other direction which is non-tangential to the criti- 
cal line, then b 2 <C a and c™ 9 is inversely proportional 
to the square root of the distance to the TCP [ToL |48| . 
This demonstrates that the critical region being defined 
by the properties of the c n for n > 1 is elongated along 
the 0(4) critical line. 

The thermodynamic properties and critical behavior 
near TCP arc well described by the mean-field theory up 
to logarithmic corrections, because in this case the upper 
critical dimension is three. 
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For the non-zero external field (non-zero quark mass) 
the three-critical point is turned to the CEP. Following 
the same mean-field analysis (see also Ref. |48|) we ob- 
tain the following leading singular behavior along the 
linear continuation of the phase-coexistence line to the 
crossover region and near to the CEP 



sing 



(46) 



For any other directions which arc asymptotically not 
equivalent to the previous one 



„sing 



(47) 



where the introduced variables v and u are linear combi- 
nations of (T — Tqep) and (/i — Hcep)- For scaling (|46|) 
the kurtosis diverges as R 
(|47| the kurtosis KT 9 



sing 
4,2 

2 



and for scaling 



"-4,2 



II 



When going beyond the mean-field approximation by 
including quantum and thermal fluctuations in the PQM 
model one expects that along the u = line the following 
scaling holds (see also |39j): 



-[(n-2)( 7 +/3)+ 7 ] 



and for any other direction 



„smg 



| u |-(n-3+^) 



(48) 



(49) 



Here, the critical exponents 7 and /3 correspond to the 
three-dimensional spin model belonging to the Z(2) uni- 
versality class [5(| ■ The kurtosis is divergent at the CEP, 
however the strength of the singularity depends on the 
direction. Approaching the TCP along the u = line 
results in 



IV. SUMMARY AND CONCLUSIONS 



We have formulated thermodynamics of the Polyakov 
loop extended quark-meson effective chiral model 
(PQM), including quantum fluctuations within the func- 
tional rcnormalization group method (FRG). We have 
solved the flow equation for the scale dependent thermo- 
dynamic potential at finite temperature and density in 
the presence of a background gluonic field. 

We have shown that the non-pcrturbative dynamics in- 
troduced by the FRG approach essentially modifies pre- 
dictions of the model derived under the mean-field ap- 
proximation. In particular, we have demonstrated quan- 
titative changes of the phase diagram leading to a shift 
in the position of the critical end point. 

We have focused on fluctuations of the net-quark num- 
ber density and calculated the first four moments near 
the chiral transition for different values of the chemical 
potential. We have indicated the role and importance of 
ratios of different cumulant moments to identify the de- 
confinement and chiral phase transitions. We have also 
discussed predictions of the Landau and scaling theories 
on the critical behavior of the net-quark number density 
fluctuations and their higher moments in the vicinity of 
the chiral phase transition. 

The extension of the FRG method proposed here to 
account for the coupling of fermions to the background 
gluon fields within the quark-meson model is of relevance 
to understand effectively thermodynamics of the QCD 
near the chiral phase transition. 



rising 

n 42 



I -2(7+0) 



-2-7+Q 



whereas for any other direction 



paws 
n 4,2 



(50) 



(51) 
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